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Abstract 

The differential cross section of da(e + e~ -> e + e~ D* ± X)/dP T was measured using a soft-pion 
analysis of D* ± -> tt±D°(D°) at TRISTAN. The average y/s was 58.1 GeV and the integrated 
luminosity used in this analysis was 198 pb" 1 , respectively. 



X 



1 Introduction 

Multi-hadron production in two-photon processes is qualitatively described using the vector-meson dom- 
inance model (VDM), the quark-parton model (direct process) |Q, and resolved photon processes |^|, 
|3|, [|, ^| |7j. In order to understand them quantitatively, charm production is a good proof, because 
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the theoretical calculations have been completed to higher order (order a s ) || , 9 . In addition, cut-off 
parameters, such as P™ m , are not necessary in the theoretical calculation JlOj, [llj] and the background 
from the VDM is considered to be very small. Studies of charm-meson production should provide useful 
information concerning the gluon density in a photon, the current charm quark mass, and the intrinsic 
Pt distribution of the partons in resolved photons. The statistics of the previous experiments were too 
poor to determine these parameters jl2[ [jj, H. 

In a previous measurement by TOPAZ using the decay mode D* + — ► TrfD°(D° — > K ~ir + X) as well 
as its charge-conjugation mode (CC) an excess in charm production compared with the theoretical 
expectation was shown |9f|. Super-symmetric particle production (£) is one possible interpretation of 
this experimental excess yjl §@@§. However, the accuracy of the experiment was not good 
enough to justify this assumption. For the above-mentioned reasons, more precise measurements of 
charm production in two-photon processes have been awaited. 



2 Analysis 

2.1 Soft-pion analysis of D*+ -> n+D° 

The decay Z3* ± — * nfD°(D°) is characterized by its small Q value. The maximum transverse momentum 
of (soft-pion or n s from hereafter) with respect to the D° (from hereafter, the descriptions of the charge 
states include their charge conjugations) flight direction is only 40 MeV. Its direction can be approximated 
by the jet-axis obtained by an invariant mass algorithm Therefore, from the distribution of the 

transverse momenta of soft-pions with respect to the jet-axes (-P|L- et ), the production cross section of 
D can be precisely measured. The acceptance of this method was proven to be one order higher than 
that of the exclusive reconstruction of D*^ decay We thus tried this method for two-photon events, 
and obtained the differential cross section of D*^ as a function of Pt, the transverse momentum with 
respect to the beam axis. Higher statistics and lower systematic error measurements were expected using 
this analysis. 



2.2 Data sample 

The data were obtained by the TOPAZ detector at the TRISTAN e+e" collider of KEK. The details 
concerning the TOPAZ detector can be found elsewhere J22|. The integrated luminosity of the event 
sample used in the analysis was 198 pb _1 . The average y/s was 58.1 GeV. The trigger conditions were 
as follows: more than two tracks with Pt > 0.3 ~ 0.7 GeV and an opening angle > 45 ~ 90 degrees 
(depending on the beam condition); the neutral energy deposit in the barrel calorimeter be greater than 
2 ~ 4 GeV; or that in the endcap calorimeter be greater than 10 GeV. 

The event selection criteria for the two-photon process were tighter than in the exclusive D*^ analysis 
Jl5|] , because this signal suffered from the background of the VDM as well as single- photon exchange 
hadronic events. The selections were as follows: the number of charged tracks be > 4; the total visible 
energy in the central part of the detector be between 5 and 25 GeV; the vector sum of the transverse 
momenta of the particles with respect to the beam axis be less than 7.5 GeV; the absolute value of the 
sum of charges be < 3; the absolute value of the cosine of the thrust axis in the laboratory frame be less 
than 0.9, where neutral clusters detected by the endcap calorimeters were included; and no large energy 
clusters in the barrel calorimeter (E > 0.25Ei, eam ). In addition, we divided each event into two jets with 
respect to the plane perpendicular to the thrust axis at the laboratory frame; the cosine of the angle 
between two jets was required to be greater than -0.9. These cuts were determined using a Monte-Carlo 
simulation of the direct process events so as to maximize the charm event acceptance and its purity. A 
total of 14128 events were selected. 
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2.3 Particle selection 



The charged track selections were as follows: the closest approach to the event vertex be consistent within 
the measurement error; the number of degrees of freedom (DOF) in the track-fitting be > 3; and Pt be 
> 0.1 GeV. The 7 selections were: the cluster be detected by a barrel calorimeter; the energy be > 0.2 
GeV; and the cluster position be separated from any charged-track extrapolations by more than 10 cm 
at the surface of the barrel calorimeter. 



2.4 Jet-axis determination 

The jet-axis was obtained by using an invariant mass algorithm |2l"| , in which particles were merged 
together in an iterative way if their invariant mass was less than 2 GeV. In this algorithm, a soft-pion 
candidate, whose Pt jet mu st be calculated, was removed from the sample particles, i.e., the jet-axis 
was recalculated each time when a particle was selected to be the soft-pion candidate. Then, those jet 
candidates having invariant masses greater than 0.3 GeV and transverse momenta with respect to the 
beam-axis of greater than 1 GeV were selected. The former cut were to remove jets made of single track. 
The fake tracks by the patter recognition (duplicating tracks using the same space-points) were removed 
from the charged track sample. By the above algorithm, two or three jets were usually reconstructed 
(85% of the events). 

Those cut values were determined using a Monte-Carlo simulation while assuming a direct process. 
The best angular resolution of the soft-pions with respect to the jet-axes was obtained using the above set 
of cut values. Even when the jet-mass was as low as 0.3 GeV, the D° direction was calculated with the 
acceptable resolution which was proven by a Monte-Carlo simulation. The estimated angular resolution of 
the soft-pions was approximated by e 3 - 7 ~ 4 - 3PT degrees (Pt in GeV). The transverse momentum resolutions 
with respect to the jet-axis were 70±10 MeV in the Pt region between 0.1 and 0.6 GeV, a slowly increasing 
function of the soft-pion's Pt- 

The transverse momentum of the soft-pion was calculated with respect to the closest jet candidate 

The effects by the pair conversion electron and Dalitz pairs were estimated using a Monte-Carlo 
simulation. They were prove n to be absorbed in the background function in the fitting procedure which 



is described in Section 2.6.2 



2.5 Prjet distribution of soft-pions 

The Pt range of the soft-pion candidates was set to 0.1 < Pt < 0.6 GeV in order to cover the predicted 
region by the direct process. The Pt jet distribution for the above-mentioned Pt range is plotted in 
Figure [| (a). There is a clear peak at around the zero Pt jet region. 

In order to derive the differential cross section with respect to Pt, we selected the following Pt 
ranges: (b) 0.1 < P T < 0.2, (c) 0.2 < P T < 0.3, (d) 0.3 < P T < 0.4, and (e) 0.4 < P T < 0.6 GeV. The 
corresponding distributions are shown in Figures |l] (b), (c), (d), and (e), respectively. The binning was 
selected so as to give similar peak entries. The resolution of the Pt's of the soft-pions are smaller than 
these binnings. 

2.6 Background estimation 
2.6.1 Monte-Carlo simulations 

In the soft-pion analysis, background emulation using the experimental data is difficult, in contrast 
to the case in a full-reconstruction analysis, where wrong sign combinations can be used to emulate 
the background |lq| . We therefore used a Monte-Carlo method. The 77 generation was based on an 
equivalent photon approximation. For the parton generations, we used lowest-order formulas (Born 
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Figure 1: Distributions of P| - et of soft-pion candidates: (a) for 0.1 < Pt < 0.6 GeV, (b) 0.1 < Pt < 0.2, 
(c) 0.2 < P T < 0.3, (d) 0.3 < P T < 0.4, and (e) 0.4 < P T < 0.6. The dashed lines are the result of a 
Monte-Carlo simulation of the background, where the normalizations were carried out using the entries 
of P$ Jet > 0.02 GeV. The hatched histo grams are contributions from e + e — > (7) — ► D* ± X. The solid 
curves are the best- fitted functions which are described in the text. The higher lines are signals and 
backgrounds, and the lower are backgrounds. 
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approximation). In the resolved photon processes, LAC1 |6| parametrization was used for the parton 
density in the photon. For light-quark event generation, a P™ m of 2.5 GeV was used, which was the 
best-fitted value using the general event shapes of a 77 — * (multi-hadrons) obtained by the TOPAZ 
detector pT[ . In addition, the VDM Monte-Carlo events and e + e~ — > 7 — » events by the LUND 6.3 
(parton shower option) program were added. In LUND 6.3, the branching ratios of D° were adjusted to 
the PDG values [23|. Hadronization of partons generated by the 77 events was carried out using LUND 
6.3 string fragmentation without a parton shower option 0]. A more detailed description together with 
a detector simulation can be found in Reference Jll| . 

The branching ratio of D* + — ► 7T+P/ was assumed to be 68.1%, which was obtained by CLEO J25|. 
The V/(V+P) ratio was set to 0.75 and a u:d:s:qq ratio of 1:1:0.3:0.1 was used. 

The results of this Monte-Carlo simulation are shown in Figures [j] (a) - (e) by the dashed lines. The 
events with D* ± -> n ± D°(W) were removed from these. They fit with the experimental data in the 
higher Pt regions. In the lower Pt region, however, since the agreement of the normalization scale 
was poor, the lines in Figure |l| were normalized using the entries of Pj-jet > 0.02 GeV. However, the 
background shape (i.e., in the higher Pt jet region) agreed well. The contributions from e + e~ — > 7 — > 
are shown in Figures |] (a) - (e) by hatched histograms, and are considered to be negligible. 

2.6.2 Fit of the P| Jet spectra 

There is a difference between the experimental and Monte-Carlo data in the lower Pt region, especially 
for the normalization scale. We conclude that this is due to the facts that the cc cross section in the 
Monte-Carlo program was not correct (also, P T nm for light quark events was not correct) and that the 
higher-order effects (order a s ) were not included (for both cc and light quark-pair events). The small 
differences appeared in the slope in the Pt jet spectrum. We therefore multiplied the smoothed Monte- 
Carlo background by a first-order polynomial (a+bP^ - et ). The coefficients (a and b) of the polynomial 
were set to be free parameters. The signal shapes were assumed to be exponential functions oc e~ p, rlP . 
The (3 at each binning was obtained by fitting Monte-Carlo signal events, and was used as a fixed 
parameter in the fitting procedures, because the soft-pion angular resolution seemed to be consistent 
with the prediction of the detector simulation. 

The results of the signal and background estimations are shown in Figures [j] (a) - (e) by the solid 
lines. The peak entries obtained by the above-mentioned procedures are 50 ± 32, 126 ± 26, 83 ± 20, and 
113 ± 29, in the 0.1 - 0.2, 0.2 - 0.3, 0.3 - 0.4, and 0.4 - 0.6 GeV P T regions, respectively. In total, we 
obtained 372 ± 54 events. 

The main systematic error of the analysis is due to this fitting procedure, since we did not know the 
exact shapes of the background. We tried various background shapes, such as (1 + aP T ^ et + bP T j et + 

cPhet)' 1 b y ALEPH H and 0- + i P T,jet/ a ) b y 1 b y TOPAZ @. We defined the systematic errors 
of the fitting procedures to be the differences in the obtained peak entries by the various background 
functions. In addition, we tried to simulate various sets of the Monte-Carlo events by changing the 
generation parameters described in the previous subsection, and carried out the same fitting procedures. 
By these procedures, we obtained systematic errors of 12, 23, 18, and 9% in the 0.1 - 0.2, 0.2 - 0.3, 0.3 - 
0.4, and 0.4 - 0.6 GeV Pt ranges, respectively. 



2.7 Experimental checks of the data sample 

We carried out some miscellaneous checks on these peaks. In order to check detector bias, the charge 
states and the cos^ 7rs distribution were tested. The background-subtracted P T - et distributions are shown 
in Figures |2| (a) - (d): (a) for 7r+, and (b) n~ , (c) cos# s > 0, and (d) cos# s < 0. If D^'s had originated 
from single-photon exchange events (for example cascade decays of bb events which were produced at low 
angles with respect to the beam axis), some asymmetries may appear in the charge states, or the cos8 7ts 
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Figure 2: Distributions of Pf, - et : (a) for 7r+, (b) n~ , (c) in the forward direction, (d) in the backward 
direction, (e) for positive charge asymmetry, and (f) for negative charge asymmetry. The data points are 
background-subtracted signals using the fitting procedure described in the text. The solid curves are the 
best-best fitted signal lines. 
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distribution. These are shown in Figures |^ (e) and (f): (e) charge of ir s multiplied by cos# s be > 0, and 
(f) that be < 0. The pairs (a) and (b), (c) and (d), and (e) and (f) seem to be consistent with each 
other. The numbers of signals in these plots are 190, 182, 202, 169, 168, and 203, respectively, where 
the statistical errors of them are about 34. In the case of single-photon exchange events (bb cascade), 
the asymmetry of 1:2 should be observed between Figure || (e) and (f). In addition, the same checks on 
the higher Pt events (0.3 < Pt < 0.6 GeV) were carried out. The peak entries corresponding to the 
definitions (a)-(f) were 97, 99, 89, 107, 102, and 94 with the statistical errors of about 20, respectively. 
They are also consistent with each other. We thus conclude that these events did not originate from 
single-photon processes. 

3 Production cross section 

3.1 Systematic ambiguities 

This analysis was based on momentum measurements of inclusive pions and jet-axis reconstructions. The 
former part was well defined and the systematic errors were well known from a previous physics analysis 
using the TOPAZ detector [H . The largest systematic error in the tracking algorithm appears especially 
in the 1ow-Pt region, i.e., the 0.1-0.2 GeV region. The systematic errors were estimated using hadronic 
events via single-photon exchange processes and the LUND 6.3 Monte-Carlo simulation, since this Monte- 
Carlo method was known to reproduce many experiments. The thus-obtained systematic errors due to 
tracking were 7 and 3% at the 0.1 - 0.2 and 0.2 - 0.6 GeV Pt regions, respectively. 

The systematic error due to the jet-axis determination were evaluated by changing the maximum 
value of the invariant mass and the transverse momentum cut in selecting jets by 10%. The errors were 
11, 10, 6, and 3% in the 0.1 - 0.2, 0.2 - 0.3, 0.3 - 0.4, and 0.4 - 0.6 GeV P T ranges, respectively 

The systematic errors in the event selection were estimated in the same way by changing the cut 
values, such as the visible energy and missing Pt, by 10%, and the sum of charges by one. Those due 
to the hardware triggers were estimated using a trigger simulation program. For a charged trigger, we 
added 5% of accidental hits in the tracking detectors. The acceptance was increased by 2 ± 1.4% by this. 
The percentage of neutral triggered events was only 5.6%. Therefore the systematic errors caused by the 
energy calibrations of calorimeters and the summing amplifiers noises were considered to be negligible 
small. In total, the systematic errors were estimated to be 7%. 

In total, the systematic errors in determining the production cross sections of the soft-pions were 19, 
26, 20, and 12% in the 0.1 - 0.2, 0.2 - 0.3, 0.3 - 0.4, and 0.4 - 0.6 GeV P T ranges, respectively. From now 
on, the errors include both the statistical and systematic contributions. 

3.2 Soft-pion cross section 

We obtained the Pt differential cross section for soft-pions. The Pt of a soft-pion is almost proportional 
to that of the D*^, and also to that of the charm quark. In order to define the detector's sensitive area, 
we restricted the cosine of the soft-pion emission angle with respect to the beam axis to be within ±0.77. 
Then, the acceptance differences, especially between the direct and resolved photon processes, became 
small ( +5% higher for the resolved photon process). The acceptance was obtained using a Monte-Carlo 
simulation of the direct and resolved (LAC1) photon processes and it was a increasing function of Pt- 
The thus-obtained differential cross section is shown in Figure || and is listed in Table |l|, together with 
the theoretical predictions (described later). 

3.3 D*^ production cross section 

In order to compare this data with the previous results of the D*^ production at TRISTAN we 
unfolded the soft-pion cross section to that of Z3* ± . The method was as follows: 
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Figure 3: Differential cross section of the soft-pions versus Pt- The histograms are the theoretical 
predictions: the cross-hatched area is the direct process, the singly-hatched area the resolved process 
(LAC1), and the open area t pair production. 



P T range (GcV) 


Experiment 


Direct 


Dircct+LACl 


Direct+LACl+tt 


0.1 - 0.2 


105 ±72 


81 


145 


152 


0.2 - 0.3 


42.3 ± 14.2 


18.7 


27.1 


34.6 


0.3 - 0.4 


23.4 ±7.5 


6.0 


7.5 


13.6 


0.4 - 0.6 


6.4 ± 1.9 


1.6 


1.9 


4.0 



Table 1: da/dP T {soft - pion)(\cos9\ < 0.77) (pb/GeV). 
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Figure 4: Differential cross section of D* ± versus Pt- The open circles were obtained from a previous 
experiment, the open squares are from this experiment, and the closed squares are the combined values. 
The histogram definitions are the same as those given in Figure 3. 

1. We tuned the theory by iteratively changing parameters such as the charm-quark mass (current 
mass) so as to fit the experimental cross section {da{D* ± ) / dPr) as well as possible. The obtained 
current charm mass was 1.3 GeV. The details are described in Section |} 

2. Using the above parameter, we carried out a simulation and made a conversion matrix from the 
soft-pion P T to that of I?* ± . 

3. We multiplied this matrix by the experimental cross section of the soft-pions. 

The variation in the Pt spectrum of D* in the Monte-Carlo was a source of systematic error. We 
therefore assumed various Pt spectra of D*^ by changing the fragmentation function parameters in 
LUND 6.3 and tried unfolding. The systematic error of the unfolding was obtained to be 7%. Also 
there were strong correlations between the neighboring Pt binnings. The largest one was 85% between 
the binnings 0.1-0.2 and 0.2-0.3 GeV. The magnitudes of the other binning pairs were 30 ~ 40%. The 
errors due to these correlations were taken into account. The obtained differential cross section of 
is indicated in Figure ^ by the open squares. The previous results are indicated by the open circles p5| . 
Both experiments are consistent within their errors; the overlapping statistics of the present analysis are 
obvious. The experimental averages of both measurements are also indicated in Figure ^ by the solid 
squares and are listed in Table g, along with the theoretical predictions. 

4 Discussion 

4.1 Predictions of a two-photon process 
4.1.1 Lowest-order calculation 

The experimental cross section of D*^ in the region 1.6 < Pt < 6.6 GeV and \cos6\ < 0.77 was 24.2 ±5.0 
pb (Table |). We first compare this with the lowest-order calculation. A current-charm quark mass of 
1.5 GeV and LAC1 parametrization with A-jt^- = 0.2 GeV were used. The predictions of the direct and 
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Pt range 


Full 


Soft- 


Average 


Direct 


Dircct+LACl 


Direct+LACl+tt 


(GeV) 


reconstruction 


pion 










1.6 - 2.6 


22.9 ±9.7 


10.7 ±5.3 


13.5 ±4.7 


7.1 


11.8 


12.5 


2.6 - 3.6 


4.80 ±2.56 


5.03 ± 1.88 


4.95 ± 1.51 


2.42 


3.48 


4.22 


3.6 - 4.6 


3.79 ± 1.92 


3.01 ±0.99 


3.17 ±0.88 


1.05 


1.29 


2.08 


4.6 - 6.6 


1.38 ± 1.23 


1.28 ±0.40 


1.29 ±0.38 


0.34 


0.41 


1.02 



Table 2: da / 'dP T {D* ± ){\co S e\ < 0.77) (pb/GeV). 



resolved photon processes were 8.5 and 2.5 pb (11.0 pb in total), respectively, significantly lower than the 
experimental observation. 

4.1.2 Higher-order corrections 

Corrections on the order of a s in the QCD part of the theory are available in the analytic calculations 
|^|, |), but not in the Monte-Carlo calculations. Our Monte-Carlo simulation was based on a lowest-order 
calculation (LO), which was followed by string fragmentation using LUND 6.3. Since string fragmentation 
includes a parton-shower-like effect, the next-to-leading-order effect in the Pt spectrum of D*^ is already 
counted. For single-photon exchange events, the difference in the momentum spectra of D*^ between the 
LO-matrix element (i.e. two-jet only) and the parton-shower options in Lund 6.3 is only 9% at ^fs = 10 
GeV. We therefore concluded that the systematic error due to the hadronization process is 9%. We 
then carried out a next-to-leading-order correction (NLO) (Py-independently) of the direct process. The 
factor was 1.31, as described in reference ||. For the resolved photon process, due to the presence of the 
process jq — > ccq (a part of this is absorbed in the gluon density function in the resolved photon), the 
Pt dependent factors were obtained in the following way j2^] : 

1. We derived the Pr-dependent ratios between the higher- and lowest-order calculations for both the 
direct and resolved photon processes. 

2. We then calculated the ratios between these ratios of the direct and resolved photon processes. 

3. The final ratios were normalized so as to fit the total cross section of the higher-order calculations 
for the resolved photon process. 

The obtained P^-dependent correction were written as 

0.5CLR£ + 0.54, 

where P T was the Pt of the charm quark in GeV. Here, the charm quark mass (m c ) and the renormal- 
ization scale (/i) of 1.6 GeV and \i = ^/2m c were used, respectively. This correction factor was used as 
event weight in the Monte-Carlo simulations. 

Then, the prediction of the cross section in the range (1.6 < Pr < 6.6, \cos9\ < 0.77) was 15.6 pb, 
still lower than the experimental observation. 

4.1.3 Dependence on the current charm quark mass and gluon Pt distribution 

If the current charm-quark mass is significantly lower than 1.5 GeV, the total cross section of D*^ becomes 
large. For example, the Monte-Carlo package PYTHIA 5.6 used m c = 1.35 GeV |5(J. However, the cross 
section increased only at a low-Pp region i.e., < 2.6 GeV. The PYTHIA 5.6 program also included a gluon 
Pt of 0.44 GeV (Gaussian) inside the resolved photon jet. This may increase the high Pp cross section 
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of D* ± . We thus lowered m c to 1.3 GeV and added a Gaussian fluctuation in the gluon Pt distribution 
inside a photon of 0.44 GeV. The obtained cross section was 17.4 pb. 

The differential cross sections which this model predicted are given in Figures | | and Tables § 
The predictions agree with the experimental data in the 1.6 - 3.6 GeV Pt region. However, there is a 
large excess in the higher Pt regions (> 3.6 GeV). 

When we tried other resolved photon parameterization, such as DG ||, the predicted cross sections 
were significantly lower than the experimental data, even upon changing some of the parameters described 
so far. 



4.1.4 Ambiguities in the prediction 

We first estimate the ambiguities in the NLO correction. As has been described, string fragmentation 
describes a part of the NLO corrections. Therefore, the ambiguities lie in the hard process descriptions, 
such as three jets and hard bremsstrahlung from the resolved quark line. These are estimated by using the 
single-photon exchange events generated by the LUND 6.3 Monte-Carlo program. The D*^ momentum 
spectra between the LO-matrix with string fragmentation and parton shower options have been compared 
at yfs=lQ GeV. A difference of 9% was observed; this was considered to be a systematic ambiguity. 

Second, we considered the ambiguities in the charm mass, renormalization scale (/i), and gluon intrinsic 
Pt inside a photon. The ranges that we selected were m c = 1.3 ± 0.15 GeV, m c < [i < 2m c , and 
< p-j} uon < 0.44 GeV, respectively. In total, the cross-section ambiguity was obtained to be 11% for 
the lower Pt region and 9% for the higher Pt region. 

Third, the ambiguity in the equivalent photon approximation was studied. In our event generation, 
the initial photon from beam was considered to be almost real photon with the direction parallel to 
the beam axis. We here used a matrix element calculation by Kuroda [ |3l| in the event generation of 
the direct process. The differential cross section at the higher Pt region was increased by 25% by this 
effect. We considered that this value is a systematic error due to the Q 2 dependence of 77 system in the 
direct process. In the resolved photon case, we do not know how to take care of this effect. We thus 
approximated that all photons are real photons. Wc considered this systematic error due to the motion 
of 77 system is absorbed in that of the intrinsic gluon Pt inside the photon. 

Last, the threshold effect of the cc productions are commented. The Wvis distribution, i.e., the 
invariant mass of visible hadronic system, are shown in Figure [| (a) and (b) for the lower and higher Pt 
soft-pion signals. They were dominated in the region where Wvis > 5 GeV. The real W 77 was expected 
to be distributed above this value. We therefore concluded that the threshold effect was small in our 
detection region. We did not consider about the systematic ambiguity due to the threshold enhancement 
by this reason. 

The experimental excess at the high-P^ region (Py(D* ± ) > 3.6 GeV) with respect to the theoretical 
prediction of the direct plus resolved (LAC1) photon processes becomes a 2.9er effect. 



4.2 Possibility of t pair production 

One of the exciting possibilities to explain this excess is in terms of t pair production. A light t at the 
TRISTAN energy region is still possible if the mass difference between t and 7 is small [[f6| |l7], [l8|, |l^] . In a 
previous paper j20) , in order to explain the high Pt excess of D* ± , i pair production was introduced with 
mi~lb GeV and m 7 =12.7 GeV. These predictions are shown in Figures |^ and |] by open histograms. The 
experimental cross sections agree very well with these assumptions. The value of determines the cross- 
section excess and that of rn 7 the event shapes, such as the thrust and missing Pt distributions. These 
two distributions are shown in Figures ^| (a)-(d). Here, we had problems in translating the number of 
events into cross sections. We thus plotted the observed number of events with the theoretical predictions. 
In addition, our Monte-Carlo generator did not include any higher order (a s ) effect, such as three-jet 



11 



Figure 5: Various event shape distributions for the lower (Pt < 0.3 GeV) and higher (Pt > 0.3 GeV) 
Pt events: (a) thrust distribution in the lower Pt events, (b) that in the higher Pt region, (c) missing 
Pt distribution in the lower Pt region, (d) that in the higher Pt region, (e) \cos9\ distribution of the 
soft-pions in the lower Pt region, and (f) that in the higher Pt region. The histogram definitions are the 
same as those given in Figure 3. 



12 



events. One should therefore be cautious about them. The real thrust predictions may be softer and that 
of the missing Pt distributions may be harder due to the presence of three-jet events. 

From these figures, the experimental excesses arc considered to be high-thrust events and high missing- 
Pt events. In order to minimize this discrepancy in the missing Pr-distribution, although we must set 
a larger mass difference (m^ — rrij), it may conflict with the recent search for t pair production by the 
VENUS experiment f§. 

If the excess is due to t pair production, the angular distribution with respect to the beam axis must 
be different from that of the two-photon processes. The results are also shown in Figures || (e) and 
(f). However, the t assumption does not show good agreement. The excess events dominated in the 
low-angle region. We therefore can not strongly conclude that the excess in events originated from t pair 
production. There is a possibility of a new high-Pp process in the two-photon interaction. If this is true, 
it would be a serious background in future linear-collider experiments. 

In addition, we checked three more distributions, i.e., Wvis, the rapidity of the hadronic system, and 
the azimuthal angle between the missing Pt and the soft-pions as shown in Figures ^ (a)-(f). These 
distributions does not fit to the t assumption very well. We definitely need more studies on this process. 
Especially the electron and/or K- meson inclusive studies would be more powerful because of higher 
statistics. 

4.3 Tagged events 

In a part of the data set (integrated luminosity of 90 pb _1 ), there were forward calorimeters (FCL; made 
of BGO) which covered the polar angle region at between 3.2 and 13.6 degrees We analyzed tagged 
events using the FCL with the same analysis. We observed 35.5 ± 12.7 soft-pions, where the lowest-order 
direct process predicts ~ 26.7 events. The tagged events are therefore consistent with the expectation. 

4.4 Comparison with other experimental results 

The charm measurement using the inclusive electron method was pioneered by the VENUS group j34|. 
Their observed production rate was consistent with a theoretical prediction by the direct and resolved 
(LAC1) processes. However, the statistics are especially low in the high-P^ region. The expected t signal 
rate was only a few events. We thus concluded that the VENUS measurement was not sensitive to the 
t signal, and was consistent with our observation. A high-statistics study using high Pt electrons is 
awaited. 

Recently, a t search was carried out by the VENUS group |32]| . They have shown an upper limit 
of ~ 12.7 GeV for the case when = 15 GeV. This search marginally conflicts with our t assumption. 

5 Conclusions 

We have measured the differential cross section (da(e + e~ — > e + e~ D*^ X) / 'dPr) ■ The P/* ± s were identi- 
fied by the transverse momenta of soft-pions with respect to the jet-axes. The average ^fs was 58.1 GeV 
and the integrated luminosity of the event sample was 198 pb , respectively. We obtained 372 ± 54 
D , which represents the highest statistics so far obtained. We compared the measured cross section 
with the theoretical predictions. 
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Figure 6: Various event shape distributions for the lower (Pt < 0.3 GeV) and higher (Pt > 0.3 GeV) 
Pt events: (a) Wvis distribution in the lower Pt events, (b) that in the higher Pt region, (c) rapidity 
of hadronic system in the lower Pt region, (d) that in the higher Pt region, (e) cosine of the azimuthal 
angle between the missing Pt and the soft-pion in the lower Pt region, and (f) that in the higher Pt 
region. The histogram definitions are the same as those given in Figure 3. 
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